Excitation Spectra of Structurally Dimerized and Spin-Peierls Chains in a Magnetic 

Field 



Weiqiang Yu and Stephan Haas 
Department of Physics and Astronomy, University of Southern California, Los Angeles, CA 90089-0484 

(February 1, 2008) 

The dynamical spin structure factor and the Raman response are calculated for structurally 
dimerized and spin-Peierls chains in a magnetic field, using exact diagonalization techniques. In 
both cases there is a spin liquid phase composed of interacting singlet dimers at small fields h < h cl , 
an incommensurate regime (h c i < h < h C 2) in which the modulation of the triplet excitation spectra 
adapts to the applied field, and a fully spin polarized phase above an upper critical field h C 2- For 
Oh. structurally dimerized chains, the spin gap closes in the incommensurate phase, whereas spin-Peierls 

chains remain gapped. In the spin liquid regimes, the dominant feature of the triplet spectra is a 
one-magnon bound state, separated from a continuum of states at higher energies. There are also 
■ indications of a singlet bound state above the one-magnon triplet. 
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I. INTRODUCTION 



In the absence of a magnetic field, structuralls4.dimerized spin chains, such as have a magnetic response 

similar to spin-Peierls chains, such as CuGeOsoQ. In both cases, there is a one-magnon triplet bound state at energy 
' O ■ transfers ui = A z , where A z measures the singlet-triplet spin gap, followed by a two-magnon continuum of states 
with an onset at higher frequencies, u> = 2A z .LTThe one-magnon bound state arises fronx-the confinement of soliton- 
antisoliton pairs due to an effective attractive potential caused by the lattice dimerization.Elu In the case of spin-Peierls 
chains, there is additional coupling of the spins to the elastic degrees of freedom of the lattice, causing the softening 
of a phonon mode.0cl In contrast, the atomic positions of structurally dimerized chains are completely locked, ruling 
out a feedback between the spins and the phonons. It is this magneto-elastic feedback in the spin-Peierls compounds, 
| which allows the spin gap to remain open even at large magnetic fields^ while in the structurally dimerized chains, 
■ the gap closes beyond a critical magnetic field strength, h&, due to the deconfinement of the soliton-antisoliton pairs. 
In both cases the triplet spectra become incommensurate at h > h c \. However, we will see that for the structurally 
OS ■ dimerized chains they are gapless soliton-antisoliton continua, whereas in the spin-Peierls compounds the dominant 
one-magnon bound state remains gapped, acquiring a modulation which depends on the magnitude of the applied 

os : neide 

In the adiabatic approximation (suppressing the phonon dynamics), the effective Hamiltonian of antiferromagneti- 
cally correlated spins coupled to a crystal lattice is given by, 



N N 

H = jJ2(l + 5 r )S r -S r+1 + -J2& (1) 



where J is the Heisenberg exchange constant, S r are local lattice distortions, and K is the lattice spring constant. The 
feedback of the phonons to the spin degrees of freedom is contained in the dependence of 8 r on K . In the absence of 
a magnetic field, S r = S(— l) r . The dimgdzatiori parameter d is a constant for structurally dimerized chains, whereas 
8 oc K~^l 2 for spin-Peierls systems JlliiJO The structurally dimerized chain can be viewed as a limiting case of Eq. 
(1) with a vanishing spring constant (K — 0) and a "frozen" (inelastic) lattice modulation, which does not vary 
&3 \ with an applied magnetic field. On the contrary, spin-Paierls chains in a sufficiently high magnetic field, h > h c \, 
gain elastic energy by adjusting their lattice to the field. H£3ll3ii3 This mag-JietopeJastic distortion can be rather well 
approximated by a sinusoidal form, S r — S cos (qr), with q = 7r + 27r(S' t 2 ot /A r ).li3li3'lia With these two modulations of S r , 
the magnetic field induced transition from the dimerized to the incommensurate phase is continuous for structurally 
dimerized chains, while it is first order for spin-Peierls chains, where a jump of KS 2 /4 occurs in the elastic energy.ll3 
In this work, the spin excitation spectra in a magnetic field of these two systems are contrasted. Using exact 
diagonalization techniques on finite lattices of up to N — 24 sites with periodic boundary conditions, the triplet and 
singlet excitations are calculated, allowing a direct comparison with inelastic neutron and Raman experiments. The 
strength of this method, although restricted to relatively small lattice sizes, is the accessibility of the full excitation 
spectrum for a given cluster. As our interest is primarily in the magnetic response, we concentrate on two effective, 
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purely magnetic model Hamiltonians, Hdim and H s p, derived from H (Eq. 1). Phonon contributions, other than 
entering via the parametrization of 5 r will be neglected. 

For the structurally dimerized spin- 1/2 antiferromagnetic Heisenberg chain, the model Hamiltonian H reduces to 

N 

H dim = J^(l + cS(-l) r )S r -S r+1 . (2) 

r—l 

In a quasi-one-dimensional compound, such as KJJ.11CI3, the dimerization parameter, 5 > 0, originates from the 
alternating spacing of the spin-carrying Cu 2+ ions.EZI In the limit 5 = 1, the system is an ensemble of N/2 uncoupled 
dimers with only two energy levels per dimer, and a spin gap A 2 = 2J. In the opposite limit 6 — 0, Hdim reduces 
to the isotropic antiferromagnetic Heisenberg chain which is quasi-long-range ordered, and thus belongs to a different 
universality class from the dimerized system. For sufficiently small lattice dimerizations, a regime with sizeable 
inter-dimer interactions can be identified with the scaling properties of the massive Thirring model.El 
The effective magnetic Hamiltonian for spin-Peierls compounds is given byt!J 

N 

H sP = J^(l + <Scos(qr))S r -S r +i. (3) 
r=l 

Here, the feedback due to the interactions of the spins with the lattice phonons enters through a field dependent 
modulation of the effective nearest-neighbor exchange integral, J e ff(r) = J (1 + 5 cos (qr)), where q = ir + 2tt (Sf 0t /N) . 
In the commensurate phase (h < h c i), the modulation is fixed at q = it, and H s p is identical to Hdim- However, 
in the incommensurate regime q continuously grows from tt to 2tt, mimicking the elastic distortion of the underlying 
lattice due to the dwaamical coupling with the spins. Hence, all eigenstates of this system have a spin as well as a 
phonon component. cB 

In the subsequent section, the spin excitation spectra of structurally dimerized chains (Hdim) in a magnetic field 
are discussed, followed by a section on spin-Peierls systems (H s p). We finish with some concluding remarks. As 
we are interested in capturing the generic features, and in particular the differences, of the two physical situations 
described above, no compound specific parameters, such as inter-chain or next nearest-neighbor exchange couplings, 
are considered. Rather, our focus will be on understanding the characteristic features of the phases in the most 
elementary magnetic models. 



II. STRUCTURALLY DIMERIZED CHAIN 



It is remarkable that for some quasi-one-dimensionaLcomppunds, such as (VO)2P20j and CuGeOz, it has been 
quite difficult to establish a unique microscopic model. ErE3~Eil For example, from fits to early measurements of the 
uniform susceptibility on (VO)2P207 it has been concluded that this material is either a structurally dimerized or a 
frustrated spin-1/2 Heisenberg chain with a sizeable next-nearest-neighbor exchange coupling.E30 In both cases, a 
spin gap opens up either due to a structural or to a frustration induced dimerization, and the resulting thermodynamic 
response is quite similar for the two proposed models. Only recently, it could be shown by inelastic neutron scattering 
spectroscopy that (V O)2-P207 is indeed a structurally dimerized chain, and that it is the lattice distortion rather than 
any frustration which gives rise to the observed spin gap.EJ Similarly, there are still rather different effective parameter 
sets for 5 and J2 used in the current literature on the spin-Peierls compound CuGeO^- In one case (5 = 0.03, J2 — 0.24) 
the spin gap opens because of the lattice distortionEj, while in the other case (5 = 0.014, J2 = 0.36) the frustration alone 
is large enough to cause a spin gapE3 It is thus of particular interest to examine the full spin excitation spectrum 
of these quasi-one-dimensional materials in a magnetic field, in order to pinpoint the most relevant microscopic 
interactions. 

Let us start by discussing the phase diagram of the structurally dimerized spin-1/2 Heisenberg chain in a magnetic 
field, shown in Fig. 1. (i) for \h\ < h c %, the system is in a spin-liquid phase with a singlet-triplet spin gap A z and 
h c \ = A z ; (ii) for h c i < \h\ < h C 2 it is a gapless spin-density wave with a field-dependent modulation; and (iii) for 
1^1 > h c 2 = 2 J it is fully spin-polarized in the direction of— the applied magnetic field. To determine the dependence 



of A z on the dimerization, we use Shanks' transformation!: 



on lattices of up to N=24 sites. The asymptotic form 



of the spin gap for a given dimerization obeys a finite-size scaling relation, 

A Z (N, 5) = A Z (N = 00, 5) + A(6) exp (-T(5)N), (4) 



where the constants A(5) and T(5) are obtained from Shanks' recursive equa 



find that the initially proposed dependence, A Z (N — 00, 5) oc <5 2 ^ 3 /\/| log <5| 



In accordance with RcfE3, we 
matches our extrapolation rather 
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poorly, while the form A Z (N — oo, 8) = 2S 3 ^ 4 (shown in Fig. 1) gives an excellent fit to our data over the whole range 
of parameter space, 8 e (0, 1]. One likely reason for this discrepancy is that the initial analytical prediction is valid 
only for very small values of 8, difficult to access with a_£Lnite-size scaling procedure. Furthermore, in this regime higher 
order logarithmic corrections also become important B3 Down to rather small values of 5 (5 > 0.3), the dependence 
of the spin gap is linear to leading order, A z sa (1 + 3<5)/2, indicating that the picture of weakly interacting dimers - 
a perturbation about the limit of isolated dimers (6 = 1) - is applicable in this parameter regime. However, at low 
values of the dimerization parameter (6 < 0.3) the interactions between the dimers become increasingly important, 
leading to a deviation from the linear dependence of the spin gap on 8. 
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FIG. 1. Phase diagram of the structurally dimerized spin-1/2 Heisenberg chain in a magnetic field. 8 is the dimerization 
parameter. 



The energy gaps between the groundstate and the lowest excited singlet and triplet as a function of the lattice 
dimerization are shown in Fig. 2. The magnitudes of the gaps in this figure quantitatively resemble the thermodynamic 
limit (N — > oo), as they were obtained from Shanks' transformation.Ej In the limit of vanishing dimerization, the gaps 
disappear, indicating that the groundstate is qualitatively different for the cases of vanishing and finite dimerizatio&i 
The singlet gap is always larger than the triplet gap, and their ratio is A$/A z = 2 for most of parameter spaceB3'EII 
Hftwever, at lower values of 8 this ratio becomes smaller, possibly approaching the predicted value of -s/3 as i5 — > 
0.qE3 Unfortunately, the quality of our finite size extrapolation procedure deteriorates in this limit, and no definite 
confirmation of Ag/A z = v3 can be drawn from this study, although our data is consistent with this value. It is 
clear, however, that the region of small dimerization must be governed by a field theory such as the massive Thirring 
model with a non-linear dependence of the excitation gaps on the effective mass, which in turn is proportional to S.u 




FIG. 2. Singlet and triplet spin gaps in the structurally dimerized spin-1/2 Heisenberg chain at zero magnetic field, obtained 
by Shanks' transformation on lattices of up to N=24 sites, (i) diamonds: singlet-triplet gap, (ii) triangles: singlet-singlet gap, 
(iii) circles: ratio of singlet- singlet and singlet-triplet gap. 
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In Fig. 3, the triplet excitation spectra of the structurally dimcrized spin-1/2 Heisenberg chain are shown, along 
with the corresponding dispersion relations in the insets. These spectra were calculated by an exact numerical 
diagonalization of 18-site chains with periodic boundary conditions, combined with a continued fraction expansion to 
obtain the full dynamical response functions. Let us first examine the dynamical spin structure factor, 

S"(k,u>) = J2 K«l^|0)| 2 J(o; - E n + E ), (5) 

n 

where S| = ^= exp(ikr)S* is the projection of the spin operator parallel to the applied magnetic field, \n) denotes 
an eigenstate of the Hamiltonian with energy E n , orthogonal to the groundstate with energy Eq, and magnetization 
m = (Sf ot /N) € [0, 1/2]. In the dimerized phase (Fig. 3(a)), the dominating feature in the spectrum is the one-magnon 
bound state with an onset frequency lu — A z w (1 + 3(5) /2, well separated from a continuum of states starting at twice 
this energy. Increasing the magnetic field, the one-magnon bound state moves down to lower energies, and eventually 
the gap closes at h = h c \ (Fig. 3(b)). Beyond h c \, the soliton-antisoliton confinement potential is thus overcome 
by the magnetic field, and the bound state decays into a low-energy two-spinon continuum, similar to the spin-1/2 
Heisenberg chain. In addition, there are continua of states at higher energies. In particular the lowest continuum 
(starting at u = A z ) carries most of the spectral weight. At the onset of the incommensurate phase (h = h c i, Fig. 
3(b)) the phase space of this low-energy continuum is strongly restricted, and therefore its width is small. The reason 
will become clear from the discussion below in terms of the corresponding spinless fermion picture. As the applied 
field is increased from h c \ to h C 2, the width of the low-energy continuum grows, and the wave vector of the dominant 
infrared divergence moves continuously from q — ir to q — 2ir. In Fig. 3(c), the triplet excitation spectra are shown 
at a particular magnetization, m = 4/18, corresponding to a magnetic field h ~ 1.58J. Clearly, the modulation at 
this field is incommensurate, and the magnetic unit cell is enlarged by approximately a factor of two with respect to 
its size at zero field. Furthermore, the phase space for triplet excitations is reduced with increasing magnetic field, 
leading to an overall loss of spectral weight at higher fields. Finally, close to h c2 two triplet bands emerge, split by a 
dimerization gap, A± = 26 (Fig. 3(d)). Low-energy spectral weight away from long wavelengths disappears, and the 
low- frequency dispersion approaches wafc 2 , characteristic for ferromagnetism. 
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FIG. 3. Triplet excitation spectra of the structurally dimerized spin-1/2 Heisenberg chain in a magnetic field, calculated on 
an 18-site lattice with S = 0.4. The lowest curve is the dynamical spin structure factor at momentum transfer k = 0, and the 
highest at fc = 7r. (a) Dimerized regime (h = 0), (b) onset of incommensurate phase (h — h c i), (c) center of incommensurate 
phase (h c i < h < h C 2) , (d) incommensurate phase at a high magnetic field, close to h C 2- The corresponding pole positions are 
shown in the insets. 

The Hamiltonian of the structurally dimerized spin-1/2 Heisenberg chain can be mapped onto a model of spinless 
fermions via the Jordan- Wigner transformation. Due to the lattice dimerization, the spinless fermion band is split 
into two parts which disperse according to 

w ± (k)/J = 1 ± v/<5 2 + (1 - S 2 ) cos 2 (fc), (6) 

giving rise to a dimerization gap, A± = 2(5, and to a total single particle bandwidth of 2J. While these dispersion 
relations are exact in the XY- limit of Hdim, they are only slightly renormalized in the isotropic Heisenberg limit for 
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sufficiently large dimerization values (5 > 0.3), as can be seen by comparing oj-(k) and to+(k) with the dispersions 
in the inset of Fig. 3(d). Furthermore, in the spinless fermion picture, the applied field corresponds to a chemical 
potential. At vanishing magnetic field, the chemical potential lies in the center of the gap between w_ and ui+. In 
order to excite an unbound pair of particles, a minimum energy of 2A± is needed. In addition, due to the attractive 
scattering between the spinless fermions, an exciton-type particle-hole bound state is formed with a minimum energy 
of A 2 (/i), dispersing at h = as 

uj z {k)/J = (1 + 5)- (1 - <5)cos(2fc)/2, (7) 

as observed in the inset of Fig. 3(a) B At small magnetic fields (h < h c i) this dominant one-magnon triplet mode 
carries most of the weight in the dynamical structure factor. Furthermore, there is a second gap (A2 = A z ) between 
the one-magnon bound state and a continuum of states which is a simple convolution of two magnons with a dispersion 
uj z (k). With increasing magnetic field (h — > h c \) the bound state moves down to lower energies, and eventually A z 
vanishes at h c \, whereas the onset of the continuum now occurs at = h c i) = A z (h = 0). Beyond h c \ the 

one-magnon bound state disappears and decays into a particle-hole continuum as the effective confining potential is 
overcome by the applied field. 

Using only u) z (k) and w±(k), the complete magnetic field dependence of the triplet spectra can thus be understood 
qualitatively within a simple rigid-band picture. In the incommensurate phase, the chemical potential moves into 
the lower band, uj-(k). The continuum of states at low energies arises from two-particle excitations within w_(fc), 
whereas the continua at higher frequencies stem from processes involving interband scattering. The modulation wave 
vector q, corresponding to the applied magnetic field h, is obtained from the solution of h = oj-(q). At fields slightly 
above h c \, the phase space for intraband scattering processes are restricted within the lower band which is almost full. 
This is the reason for the narrow width of the low-energy continuum in Fig. 3(b), just at h = h c i- At large magnetic 
fields, h < h C 2, the phase space of triplet excitations is exhausted, and the dynamical structure factor traces out the 
single particle bands of the spinless fermions (Fig. 3(d)). 

Because the spectra in Fig. 3 were obtained on finite-size lattices, there are no true branch cuts in S zz (k,to). 
Rather, discrete sets of poles appear where continua are expected to emerge in the thermodynamic limit. In order 
to distinguish bound state poles from sets of poles which become part of a continuum in the thermodynamic limit, a 
finite-size scaling analysis of the individual pole positions and weights is necessary. From such an extrapolation, using 
chains with N = 4, ... , 24 sites, we find that the bound state with oj z (k) (Fig. 3(a)) is indeed stable, whereas the other 
poles in the spectrum merge into continua as the lattice size is increased to infinity. At smaller values of the lattice 
dimerization (S < 0.3), the general features of the triplet spectra are the same as discussed above. However, once 
the bandwidths of u>±(k) become larger than the dimerization gap separating the two bands, A±, the higher-energy 
continua, which are separated from each other for larger dimerizations, begin to overlap. Using the exact dispersion 
expressions for the XY limit (Eq. 6), this crossover occurs at S c = 1/3 w 0.3, consistent with the deviations from the 
weakly interacting dimer picture observed in our numerical data. 




FIG. 4. Raman spectra for a 20-site spin-1/2 structurally dimerized Heisenberg chain in a magnetic field with (b) S — 0.1 
and (c) 5 — 0.4. The integrated weight for these two cases is shown in (a). The spectra are for magnetizations m — (lowest 
curve) up to m = 9/20 (top curve). 

Let us now turn to the spin excitation spectra of the structurally dimerized spin.ehain in a magnetic field, as they 
are probed by Raman scattering measurements. Within the Loudon-Fleury theoryE2l - assuming resonant scattering 
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- the effective A\ g Raman operator of a one-dimensional spin system is proportional to ^ r S r ■ S r +i. Setting the 
proportionality constant equal to one, the dynamical Raman response function takes the form 

J M = X! Sr ' S r+i\0)\ 2 S(uj - E n + E ). (8) 

n r 

In the following discussion of I(u>), terms in the Raman operator which are proportional to the Hamiltonian are 
omittedo In Fig. 4 the Raman spectra for 20-site chains are shown as a function of the magnetic field. At zero 
magnetization (h < h c i), one singletJhound state is expected with a gap A s , followed by a continuum of excitations at 
higher energies, involving 4 spinons.o In the regime of large dimerization (8 > 6 C ), the singlet gap is at A s — 2A Z ~ 
(1 + 3<5) J (Fig. 4(c)). This bound state is best understood by considering the limit of complete dimerization (<5 = 1). 
The corresponding groundstate at h — is a product of singlet dimers with energy — 3J/2 per dimerEj Two such 
dimer singlets can be excited by the Raman operator into a 4-site singlet state. The energy difference, A s , between 
these states is approximately equal to (1 + 36) J as long as the dimer-dimer interactions are sufficiently small (S > 5 C ). 
In the limit of complete dimerization, one obtains exactly A s = 4J. In the opposite limit, the singlet bound state 
moves to lower energies (Fig. 4(b)), and most of the spectral weight is transferred into the zero-frequency peak, not 
shown here. Apart from the bound state at A s , there is a continuum of excitations at higher energies. Evidently, 
the spectra in Fig. 4 are plagued by severe finite-size effects, such that it is particularly difficult to distinguish by 
inspection the precursors of continua from emerging isolated bound states.EJ However, from the scaling behavior of the 
individual pole positions and weights we conclude that our data at zero magnetic field is consistent with the picture 
of an isolated bound state at A s , followed by a continuum of states above a threshold u> > A z in the thermodynamic 
limit. Beyond h c i, the bound state disappears, and the onset frequency of the continuum increases from A z at h c \ 
up to 2 J close to h c2 . In Fig. 4(a) the integrated weight of the Raman spectrum, W = J <1loI{lu) : is plotted as a 
function of the magnetization. W becomes very small as S — > 0. Furthermore, with increasing magnetic field, W(m) 
decreases more rapidly - with a purely concave shape - in the case of small dimerization (Fig. 4(b)). However, this 
subtle difference is most likely of little experimental relevance because important compound specific contributions, 
such as frustrating longer range interactions or interchain coupling, have been neglected in this discussion. 



III. SPIN-PEIERLS CHAIN 




FIG. 5. Magnetization curves of the spin-1/2 structurally dimerized Heisenberg chain (circles) and of the spin-Peierls chain 
(triangles) with 5 — 0.4. At h c i, the magnetization of the spin-Peierls system jumps discontinuously from zero to a finite value, 
indicating a first order transition, whereas in the structurally dimerized chain the transition into the incommensurate regime 
is continuous. 

In contrast to the structurally dimerized chains, spin-Peierls chains have a gapped incommensurate phase as their 
lattice distortion adapts magneto-elastically to the applied field. For a self-consistent treatment of this incommensurate 
regime, the phonon and the spin degrees of freedom have thus to be treated on an equal footing. In large scale numerical 
studies of the full adiabatic spin-phonon Hamiltonian H (Eq. 1), it has been shown that the structural distortion 
of the lattice and the modulation of the local magnetization have the shape of solitons, natural for one-dimensional 
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systems Jla€3 As the focus of this w.ojk is the magnetic response, an effective parametrization for the lattice distortion 
of the form S r = 5 cos(qr) is usedJ13 instead of treating the elastic part of H self-consistently. Especially at higher 
magnetic fields, this form gives results almost identical to the self-consistent treatment of the lattice dynamics for. 
observables such as the local magnetization. Furthermore, we have verified that other (solitonic) parametrizationsE3 
yield spin excitation spectra which are only minutely different to those with a sinusoidal lattice distortion. 

In Fig. 5, the magnetization curves, m(h), of the spin-1/2 structurally dimerized Heisenberg chain and of the spin- 
Peierls system are shown. They were calculated by numerical diagonalization of 22-site chains. The magnetization 
of the structurally dimerized chain has a plateau between h — and h c i — A z . m(h) then rises continuously from 
zero at h c \ up to rn = 1/2 at h C 2 = 2 J. The particle-hole symmetry of the corresponding spinless fermion band is 
reflected by the point symmetry of m{h) about its midpoint. In contrast, the magnetization of the spin-Peierls chain 
jumps discontinuously at h c i from zero to a finite value. The position of h c \ and the magnitude of the discontinuity 
depend on the lattice spring constant K. Therefore, a precise determination of these quantities is beyond the realm of 
our effective - purely magnetic - theory using H s p. However, from self-consistent calculations it has been concluded 
that the second order transition line, h c i(S), for the structurally dimerized chain (Fig. 1) becamfts-jfirst order in the 
spin-Pcicrls case, and moves down towards lower fields with increasing lattice spring constant EOEj 
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FIG. 6. Triplet excitation spectra of the spin-1/2 spin-Peierls chain in a magnetic field, calculated on an 18-site lattice with 
S — 0.4. The lowest curve is the dynamical spin structure factor at momentum transfer k = 0, and the highest at k = it. (a) 
Dimerized regime (h = 0), (b) onset of incommensurate phase (h = /i c i), (c) center of incommensurate phase (h c \ < h < h C 2) 
, (d) incommensurate phase at a high magnetic field, close to h C 2- The corresponding pole positions are shown in the insets. 

Let us now turn to the triplet excitation spectrum of the spin-Peierls chain in a magnetic field, shown in Fig. 6. For 
h < h c i, the spin response is identical to that of the structurally dimerized chain (Fig. 6(a)). The lattice dimerization 
gives rise to a scattering potential peaked at q = 2kp = n. This attraction between pairs of particles at opposite ends 
of the Brillouin zone leads to the formation of a one-magnon bound state. In the incommensurate phase (Fig. 6(b-d)), 
the lattice distortion, q = n + 2nm(h), adapts to the magnetic field, thus supporting a bound state with a minimum 
energy u> — A z (h) ^ A z (h — 0) and a modulation q due to an effective field dependent potential, V e ff(q) ~ 8cos(qr). 
At higher energies, there are continua of states. In contrast to the dimerized phase, in the incommensurate regime the 
second gap, between the onset of the bound state and the onset of the lowest continuum, is smaller than A z (h). This 
can be understood within the spinless fermion picture: as the magnetic field is increased beyond h c i, the corresponding 
Fermi wave vector moves away from 7r/2. Opposite to the case of the structurally dimerized chain, the scattering 
potential adapts to the changing magnetic field, such that there is always an instability at the Fermi level due to 
particle- hole scattering with momentum transfer q = 2kp(h) = n + 2nm(h). Furthermore, in the incommensurate 
phase the chemical potential is offset from the center of the gap. The energy difference between [i and the lower edge 
of the gap is A_(ft,), and between n and the upper edge it is A + (h), shawn in Fig. 7, where the values of A_ and 
A + have been evaluated by acting with S~ and S + on the groundstate.tL2l The onset of the one-magnon bound state 
occurs at a higher energy, A z (h) > (A^(h) + A + (h))/2. Only in the dimerized phase, it is found that A z = A_ = A + , 
and thus A z = (A_ + A+)/2 because of the additional particle-hole symmetry for the special case of half-filling. The 
magnetic field dependence of the triplet spectra in Fig. 6 follows exactly this picture. For example, at magnetization 
to = 4/18 (Fig. 6(c)), the onset of the one-magnon bound state is at u> — 0.82 J, whereas the lowest continuum of 
states starts at A_ + A + = 0.39J + 0.69 J = 1.08 J, thus separating the bound state from the continuum by a second 
gap of A 2 = 0.26 J < A z {h) = 0.82 J. 



7 



Similar to the structurally dimerized chain, the phase space of triplet excitations is gradually reduced in the 
incommensurate phase of the spin-Peierls system, leading to a loss of spectral weight in S zz (k,ui) with increasing 
magnetic field. Furthermore, the widths of the triplet bands shrink at higher fields, indicating an effective localization 
of triplets. For example, the band width W of the one-magnon bound state is drastically reduced (see inset of Fig. 
7), going practically to zero beyond m c rs 0.3. Its dispersion, u> z (k), oscillates rapidly as the modulation vector q 
approaches 2tt. The reason for this behavior is that the real-space magnetic unit cell grows with increasing magnetic 
field. At sufficiently high fields, it grows beyond the size of any finite cluster. For the parameter choice and lattice size 
we use, this happens at approximately m c . Beyond m c , the corresponding modulation of the effective nearest-neighbor 
exchange integral, J e //(r) = J(l + Scos(qr)), has only one minimum in the finite chain where triplets are trapped, 
leading to a "smearing" of u>(k) in momentum space (Fig. 6(d)). While this localization of magnons is obviously an 
artefact of the finite cluster calculation, it may be realized in mesoscopic chains, as soon as the size of the magnetic 
unit cell exceeds the mesoscopic length scale. Also, such a localization can easily be stabilized by a pinning of the 
distortion, to which a physical system is highly susceptible as the modulation grows toward infinity. In this case, a 
real-space picture of (almost) localized triplets is most appropriate at high magnetic fields, close to h c2 . 
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FIG. 7. Singlet and triplet gaps in the incommensurate phase of the spin-Peierls chain with 5 = 0.4, calculated on a 20-site 
lattice, (i) circles: triplet gap parallel to the applied field, (ii) squares and diamonds: triplet gaps perpendicular to the applied 
field, (iii) triangles: singlet gap. 

The singlet gap A s , shown in Fig. 7, corresponds to the onset of the finite frequency Raman spectrum (Fig. 8). It 
is always larger than the triplet gap, and has a similar dependence on the magnetization. From an examination of 
the finite-size scaling behavintjof the poles in the Raman spectrum, there appears to be a low-energy singlet bound 
state for all fields h < h^u'cB In particular, an analysis of the 5 — 0.4 Raman spectra (Fig. 8(c)) suggests that 
in the thermodynamic limit the two lowest poles merge into one, and their spectral weight extrapolates to a finite 
value, thus indicating the existence of a two-magnon bound state. For smaller dimerizations (such as S = 0.1 in 
Fig. 8(b)) it is difficult to determine from our finite-size data whether there is a bound state. This is consistent 
with a recent numerical study of the Raman spectrum in CuGeO^ which has considered an even smaller dimerization 
constant S ~ 0.03.E3 Here it was argued that the strong magnon-magnon interactions destabilize the singlet bound 
state. However, as seen in Figs. 8(b) and (c), the Raman excitation spectra are qualitatively rather similar for these 
two parameter choices. Consider for example the spectra at the onset of the incommensurate phase (to = 1/20), 
shown in the second lowest curves of Figs. 8(b) and (c). At low energies, there is a pair of poles (the second pole is 
not visible for S = 0.1), separated from a set of poles at higher energies with a close spacing. As discussed above, the 
low- frequency poles merge in the thermodynamic limit, whereas the high-frequency set of poles appears to evolve into 
a continuum of states. We therefore suspect that singlet bound states may exist at the lower edge of the spectrum 
for any finite 5, but much larger clusters may be needed for a numerical confirmation of the singlet bound state at 
small values of 5. Furthermore, the total finite- frequency spectral weight W(m), shown in Fig. 8(a), has the same 
dependence on the magnetization for both choices of <5, indicating that there should be only one massive field theory 
applicable for the whole range 5 € (0, 1]. 
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FIG. 8. Raman spectra for a 20-site spin-Peierls chain in a magnetic field with (b) 5 = 0.1 and (c) S = 0.4. The integrated 
weight for these two cases is shown in (a). The spectra are for magnetizations m = (lowest curve) up to m = 9/20 (top 
curve) . 



IV. CONCLUSIONS 



In summary, wc have analyzed the spin excitation spectra for two distinct models in a magnetic field: the structurally 
dimcrized chain and the spin-Peierls chain. Below a critical field, ft, cl , both systems are in a spin liquid phase, composed 
of interacting singlet dimers. Above h c \, the spin gap of the structurally dimerized chain closes, whereas the spin- 
Peierls system supports a singlet-triplet gap up to h C 2- Therefore, the excitation spectra of these two models are quite 
different in their incommensurate phases. In the structurally dimerized chain, a soliton-antisoliton continuum appears 
at low energies, separated by a dimerization gap A± = 25 from a second continuum at higher frequencies. In the 
spin-Peierls chain there is a triplet bound state with an onset at A z , and a higher-energy continuum of states, starting 
at A_ + A + < 2A Z . Parts of the triplet band may thus overlap with the continuum. Common features in the spin 
excitation spectra of these two systems are (i) an incommensurate, field-dependent modulation q = ir + 2irm(h) for 
h c \ < h < h C 2, (ii) a loss of overall spectral weight with increasing magnetic field, and (iii) and full spin polarization 
beyond h C 2- 

Furthermore, there are qualitative changes in the spectra of the structurally dimerized chain, depending on the 
magnitude of the dimerization parameter 6. The region of large dimerization in the h — S phase diagram can be 
understood within the valence bond picture of weakly interacting dimers, whereas for small values of S the interactions 
between the dimers are important, reducing the ratio of the singlet to the triplet gap and increasing the bandwidths 
of the spectral features. 

The triplet excitation spectra of the spin-Peierls chain in the incommensurate phase contain a one-magnon bound 
state with an onset at A z , and a soliton-antisoliton continuum with an onset at higher frequencies, ui = A_ + A + . 
The actual dimerization strengths of real spin-Peierls compounds are typically smaller than the values we have studied 
here, which were chosen to improve numerical stability. However, the spin excitation spectra obtained for H s p arc 
qualitatively similar for the whole range of <5 we were able to study. 

We wish to thank A. Honcckcr, B. Normand, G. Uhrig, and S. Wessel for useful discussions, and acknowledge the 
Zumberge Foundation for financial support. 
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